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Abstract: Changes in atmospheric aerosols have profound effects on ecosystem productivity, vegetation 
growth and activity by directly and indirectly influencing climate and environment conditions. However, 
few studies have focused on the effects of atmospheric aerosols on vegetation growth and activity in the 
vulnerable arid and semi-arid regions, which are also the source areas of aerosols. Using the datasets of 
aerosol optical depth (AOD), normalized difference vegetation index (NDVI) and multiple climatic 
vatiables including photosynthetically active radiation (PAR), surface solar radiation (SSR), surface air 
temperature (TEM) and total precipitation (PRE), we analyzed the potential responses of vegetation 
activity to atmospheric aerosols and their associated climatic factors in arid and semi-arid regions of Asia 
from 2005 to 2015. Our results suggested that areas with decreasing growing-season NDVI were mainly 
observed in regions with relatively sparse vegetation coverage, while AOD tended to increase as NDVI 
decreased in these regions. Upon further analysis, we found that aerosols might exert a negative influence 
on vegetation activity by reducing SSR, PAR and TEM, as well as suppressing PRE in most arid and 
semi-arid regions of Asia. Moreover, the responses of atmospheric aerosols on vegetation activity varied 
among different growing stages. At the early growing stage, higher concentration of aerosol was 
accompanied with suppressed vegetation growth by enhancing cooling effects and reducing SSR and PAR. 
At the middle growing stage, aerosols tended to alter microphysical properties of clouds with suppressed 
PRE, thereby restricting vegetation growth. At the late growing stage, aerosols exerted significantly 
positive influences on vegetation activity by increasing SSR, PAR and TEM in regions with high 
anthropogenic aerosols. Overall, at different growing stages, aerosols could influence vegetation activity by 
changing different climatic factors including SSR, PAR, TEM and PRE in arid and semi-arid regions of 
Asia. This study not only clarifies the impacts of aerosols on vegetation activity in source areas, but also 
explains the roles of aerosols in climate. 
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1 Introduction 


Drylands, including arid and semi-arid regions, cover about 41% of the world's land surface area 
(Huang et al., 2017). Over the past centuries, drylands have dramatically transformed in response 
to frequent human activities and climate changes (IPCC, 2013). Vegetation growth is vulnerable 
to harsh environments in drylands (Peng et al., 2011; Pauli et al., 2012; Zhang et al., 2014; Wang 
et al., 2017), and climatic factors including solar radiation, temperature and precipitation may 
have profound impacts on the vegetation dynamics in drylands (Nemani et al., 2003; Park and 
Sohn, 2010; Mohammat et al., 2013; Kim et al., 2014). 

Drylands are also regarded as the potential sources of aerosols (Shao et al., 2011). Since the 
2000s, atmospheric aerosols have increased rapidly, which would have remarkable influences on 
the surface energy budget of the Earth (Yoon et al., 2014; Hsu et al., 2017). On the one hand, 
aerosols could modify the shortwave radiation that the Earth receives by scattering and absorbing 
solar radiation (Chameides et al., 1999; Steiner et al., 2005; Wang et al., 2018). On the other hand, 
aerosols could interact with cloud droplets by serving as cloud condensation nuclei (CCN) or ice 
nuclei (IN), which would change the physical and radiative properties of cloud, and eventually 
alter the Earth's radiation balance (Chameides et al., 1999; Steiner et al., 2005; Wang et al., 2018). 
Not only that, interactions between aerosols-induced radiation changes and microphysical 
processes of clouds could further produce complex effects on climatic factors (Gu et al., 2016; 
Ezhova et al., 2018). For instance, aerosols have been reported to cool the Earth surface with a 
particular reduction in night-time temperature (Nair et al., 2011). Moreover, since aerosols are 
involved in the microphysical processes of clouds (Rosenfeld et al., 2014; Dong et al., 2019), they 
would substantially affect regional precipitation, but their effects on precipitation differ among 
situations with various types of aerosols and cloud (Hui et al., 2008; Min et al., 2009; Yan et al., 
2014). Among all the types of aerosols, dust aerosols are found to change the number of cloud 
droplets and cloud water path, thereby affecting low cloud cover, high cloud cover and water 
vapor contents, which ultimately alter the precipitation rate (Huang et al., 2006). In addition, dust 
aerosols may reduce the cloud particle sizes and even 'burn' the clouds (i.e., the heating effect; 
Yin and Chen, 2007), thus suppressing warm rain formation (Rosenfeld et al., 2001; Huang et al., 
2009; Qian et al., 2009). 

Radiative effects of aerosols could indirectly influence ecosystem productivity and CO2 
exchange through multiple mechanisms (Tie et al., 2016; Yue et al., 2017; Wang et al., 2018), 
which are found to vary among different ecosystems (Yue et al., 2017). Generally, the response of 
net ecosystem CO2 exchange to aerosol loading is negative in grasslands (including semi-arid 
grasslands), but positive in forests and croplands, possibly due to their different canopy structures 
(Roderick et al., 2001; Niyogi et al., 2004, Xie et al., 2010). 

Arid and semi-arid regions of Asia is considered as one of the major source areas of dust 
aerosols, which has emitted a total of >30% of global dust aerosols (Zhang et al., 1997; Ginoux et 
al., 2012). This region covers an area of 1.3x10’ km? stretching from central and western to 
northeastern Asia, and receives an annual mean precipitation of <500 mm (Fig. 1), which makes it 
extremely sensitive to climatic and environmental changes. For instance, Gessner et al. (2013) 
found that changes in the seasonal precipitation of arid and semi-arid regions in Asia, especially 
in Central Asia with an annual precipitation of 100—400 mm, were closely related to the evolution 
of vegetation. In addition, Zhang et al. (2016) found that increases in winter precipitation and 
summer temperature may be beneficial for vegetation greening in eastern areas of Asia. 
Mohammat et al. (2013) found that drought and spring cooling may reduce vegetation growth in 
the extremely arid regions. However, potential effects of aerosols on vegetation activity in arid 
and semi-arid regions of Asia are still poorly understood. 

In order to understand the direct effect of aerosols on climatic factors and their indirect effects 
on vegetation activity in arid and semi-arid regions of Asia, we used datasets of aerosol optical 
depth (AOD), normalized difference vegetation index (NDVI), and major climatic factors 
including photosynthetically active radiation (PAR) (surface visible solar radiation with a spectral 
range of 400—700 nm), surface solar radiation (SSR), surface air temperature (TEM) and total 
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precipitation (PRE) to analyze the spatiotemporal variations in the correlations between 
vegetation activity and major climatic factors that are potentially influenced by aerosols. In this 
study, we aimed not only to identify the ecological effects of aerosols in their source areas, but 
also to clarify the roles of aerosols in climate and ecosystem. 
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Fig. 1 Arid and semi-arid regions and land uses in Asia. Land uses are referenced from NASA (National 
Aeronautics and Space Administration; https://lpdaac.usgs.gov). AERONET, aerosol robotic network. 1, 
Dushanbe; 2, Issyk-Kul; 3, Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL); 4, 
Dalanzadgad. 


2 Data sources and methods 


2.1 Data sources 


NDVI data from 2005 to 2015 were obtained from MODIS (moderate-resolution imaging 
spectroradiometer)/Terra collection 5.0 Vegetation Index Products (MOD13C1-NDVI), which 
were downloaded from LAADS (level 1 and atmosphere archive and distribution system) web 
(https://ladsweb.nascom.nasa.gov/). This gridded dataset (with a temporal resolution of 16-d and 
a spatial resolution of 0.05°x0.05°) was predominantly designed to provide a globally consistent 
assessment of vegetation condition (Running et al., 1994; Justice et al., 1998; Salano et al., 2010). 
To improve the accuracy of the analysis results, we defined those pixels with a mean 
growing-season (81-304 d) NDVI of >0.1 as vegetated areas (Zhou et al., 2001). 

AOD has been widely used to develop global aerosol pattern as it comprehensively describes 
the optical properties of aerosols and quantitatively describes aerosol radiative forcing (Levy et 
al., 2013; Yu et al., 2013). In this study, we used the daily AOD data from 2005 to 2015, with a 
spatial resolution of 10 kmx10 km, from the MODIS Aqua collection 6.0 Aerosol Product 
(MYD04). In the datasets of MYD04 products, the scientific dataset 
"‘Dark_Target_Deep_Blue_Optical_550_Combined' was used for further analysis, which could 
select one type of retrievals (if any) or an average of the two retrievals that applied the dark target 
or deep blue algorithm. The ground-observation AOD was obtained from aerosol robotic network 
(AERONET) of NASA (National Aeronautics and Space Administration). CE318 sunphotometer 
was deployed in the AERONET for detecting the bands at wavelengths of 380, 440, 500, 675, 870 
and 1020 nm. AERONET measurement was generally used to validate the satellite-retrieved AOD 
since it contained the spectral data on solar radiation, from which various aerosol optical 
properties could be retrieved using an inversion algorithm (Dubovik et al., 2000). The four 
AERONET measurement sites of AOD data in arid and semi-arid regions of Asia (including 
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Dushanbe (38.55°N, 68.86°E), Issyk-Kul (42.62°N, 76.98°E), Semi-Arid Climate and 
Environment Observatory of Lanzhou University (SACOL) (35.95°N, 104.14°E), and 
Dalanzadgad (43.56°N, 104.42°E)) were used to validate the MYD04-AOD data from 2005 to 
2015 (Fig. 1). 

Since the atmospheric aerosols have affected vegetation growth primarily by altering the 
quantities and qualities of PAR or by changing climatic factors (surface temperature and 
precipitation) (Mike, 2002), the relationship between NDVI and parameters such as TEM, PRE, 
SSR and PAR were analyzed in this study. All the climatic factor data were obtained from the 
ERA-Interim reanalysis from the European Center for Medium-Range Weather Forecasts 
(ECMWF; Simmons et al., 2006). All the variables (TEM, PRE, SSR and PAR) had a spatial 
resolution of 0.125°x0.125° and a temporal resolution of 12 h (except for TEM with a 3-h 
resolution). 


2.2 Methods 

2.2.1 AOD at 550 nm wavelength 

To compare satellite-retrieved AOD with that from ground observation, it is necessary to obtain 
AOD at the same wavelength by interpolating AOD from ground observation. In this study, AOD 
and Angström exponents of 440 and 675 nm channels in the level 2.0 products were interpolated 
to a wavelength of 550 nm after cloud screening and quality control with the AERONET. The 
interpolation method can be expressed as: 


ah) Sle (1) 


where 7z(A) is the AOD of the 4 channel; P is the Angstrém turbidity coefficient; and +,,, is the 


Angstr6m exponent between wavelengths of a and b nm (here the wavelength range is 440-675 
nm). 

For the measurements at all sites collected from 2005 to 2015, we extracted the AERONET 
data at 13:30 (LST) with a 3-h time window (i.e., between 12:00 and 15:00), which just 
corresponded to the passing time of the MODIS data (Ginoux et al., 2012). AOD (550 nm) values 
from the AERONET network were then averaged during the satellite overpass time, which could 
be compared with the regional mean of MODIS retrievals within a 0.5° radius around each 
ground station (Mi et al., 2007). 

2.2.2 Linear trend in NDVI and AOD at growing stages 


AOD and NDVI at growing stages were obtained by averaging the values from 81 to 304 d every 
year. Then, linear trends in NDVI and AOD at growing stages were calculated in each pixel using 
ordinary least square (OLS) regression method: 
y=at Pxt+e, (2) 

where y is the annual mean NDVI or AOD; a is the intercept; J is the slope; x is the monitored 
year; and e is the error term. 

The slope was interpreted as a measure of linear trends in NDVI and AOD over time, and t test 
was used to determine the significance of linear trends in NDVI and AOD for each pixel. 
2.2.3 Partial correlation between NDVI and climatic factors 
To examine the potential impact of atmospheric aerosols on vegetation activity, we estimated the 
partial correlation between NDVI and climatic factors. Partial correlation is used to measure the 
linear dependence between two variables when the influence of possible controlling variables is 
excluded (Engström and Ekman, 2010). In the case of k variables (xı, X2, ...., Xx), the partial 
correlation is defined as: 
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In this study, all the datasets (AOD, NDVI and climatic factors) were preprocessed to have the 
same resolution (temporal resolution of 16-d and spatial resolution of 0.125°x0.125°). Then, we 
calculated the partial correlation between NDVI and climatic factors in two different scenarios: (1) 
the scenario a was used to calculate the partial correlation between NDVI and each climatic factor 
with the other three factors as controlling variables; and (2) the scenario b was used to get the 
partial correlation between NDVI and each climatic factor with both the other climatic factors and 
AOD as controlling variables. Significance of partial correlation could be given using a two-tailed 
t test. Finally, we compared the spatial patterns of inter-annual partial correlation between NDVI 
and four climatic factors for each 16-d period in different scenarios. 

To further explore the influence of aerosols on the correlation between NDVI and climatic 
factors at different growing stages (including early (81-144 d), middle (145-240 d) and late 
growing stage (241-304 d)), we selected five typical regions with significant NDVI change trends 
at growing stages, including the Caspian Coastal Plain (region A1), the Indus River Plain (region 
A2), eastem Mongolian Plateau (region A3), the Ordos Plateau (region A4) and eastern 
Kazakhskiy Melkosopochnik (region A5) (Fig. 1). In addition, percentages of entire region and 
five typical regions showing significantly positive correlations were also calculated for each 16-d 
period to represent the evolutions of partial correlation at different growing stages. 


3 Results 


3.1 Comparison of AERONET AOD and MODIS AOD 


In arid and semi-arid regions of Asia, comparison of AOD values from AERONET observations 
and MODIS retrievals showed general agreements at four sites (P<0.05; Fig. 2), especially for the 
Dushanbe site, which had a correlation coefficient as high as 0.79. Compared with in situ 
observation, MODIS AOD values at SACOL, Dushanbe and Dalanzadgad sites were 
underestimated, but it was slightly overestimated at Issyk-Kul site. Although there were some 
MODIS-retrieved AOD values that lay beyond the error range of AERONET observation, the 
MODIS AOD values could still represent the regional AOD characteristics in arid and semi-arid 
regions of Asia. 
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Fig. 2 Scatter plot of MODIS (moderate-resolution imaging spectroradiometer) retrieved AOD (aerosol optical 
depth) against AERONET (aerosol robotic network) AOD 
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3.2 Temporal variations of NDVI and AOD at growing stages 


NDVI climatology at growing stages from 2005 to 2015 showed that across arid and semi-arid 
regions of Asia, higher NDVI (>0.4) at growing stages was primarily found in eastern Inner 
Mongolian Plateau, northern Mongolian Plateau and eastern Kazakhskiy Melkosopochnik (Fig. 
3a), which were mainly covered by forests, croplands and grasslands (Fig. 1), whereas NDVI at 
growing stages ranging between 0.1 and 0.2 mainly scattered in the Caspian Sea Plain, Turan and 
Indus River Plains with a relatively sparse vegetation. From 2005 to 2015, about 13.2% of arid and 
semi-arid regions of Asia showed significantly NDVI increasing trends at growing stages (Fig. 
3b), which were mainly distributed in northeastern Mongolia and the Ordos Plateau, with a 
maximum increasing rate of >0.02/a. By comparison, only 4.5% of arid and semi-arid regions of 
Asia showed significantly NDVI decreasing trends at growing stages (P<0.05), with a decreasing 
rate of >0.01/a. 
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Fig. 3 Spatial pattern (a) and linear trend (b) of NDVI (normalized difference vegetation index) at growing 
stages from 2005 to 2015 in arid and semi-arid regions of Asia. Areas passing the significance test (P<0.05) are 
shown in a gray square grid in Figure 3b. 


Average AOD at growing stages from 2005 to 2015 showed that in arid and semi-arid regions 
of Asia, region with a higher AOD was mainly located in the Indus River Plain, Tarim Basin and 
the Caspian Coastal Plain, whereas a lower AOD was mainly found in northern Mongolian 
Plateau and eastern Kazakhskiy Melkosopochnik (Fig. 4a). In arid and semi-arid regions of Asia, 
regions with a higher AOD basically corresponded to those areas with a lower NDVI, except for 
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few regions such as northern Mongolian and northeastern Inner Mongolian Plateau. From 2005 to 
2015, about 50.0% of the entire study area showed increases in AOD (Fig. 4b), but only 6.8% of 
those regions were significant (P<0.05) and were mainly located in lower-elevation regions such 
as the Caspian Coastal Plain (with an increasing rate of >0.05/a). On the other hand, there were 
7.0% of total areas experiencing significantly (P<0.05) AOD decreasing trends at growing stages, 
which were mainly located at higher-elevation regions, particularly in parts of eastern Mongolian 
Plateau, with a maximum decreasing rate of >0.05/a. 
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Fig. 4 Spatial pattern (a) and linear trend (b) of AOD (aerosol optical depth) at growing stages from 2005 to 
2015 in arid and semi-arid regions of Asia. Areas passing the significance test (P<0.05) are shown in a gray 
square grid in Figure 4b. 


3.3 Linkage between NDVI and climatic factors under the influence of aerosols 


3.3.1 Effect of PAR on vegetation activity 


Partial correlation between PAR and NDVI in scenario b (i.e., taking AOD, PRE, TEM and SSR 
as the controlling variables) showed that when excluding the influence of AOD in addition to the 
other three climatic factors, area percentages with significantly positive correlation were much 
larger than those in scenario a (i.e., merely excluding the influence of PRE, TEM and SSR). This 
phenomenon was more obvious at the early growing stage than in any other stages (Table 1). 
Moreover, five typical regions were selected to show the detailed differences in partial correlation 
between scenarios a and b. For instance, at both early and late growing stages, area percentage 
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Table 1 Area percentage (%) with significantly positive correlations from two scenarios of partial correlation 
between PAR (photosynthetically active radiation) and NDVI (normalized difference vegetation index) in 
different regions at growing stages 


Early growing stage Middle growing stage Late growing stage 
Region “Scenario Scenario Scenario Scenario Scenario Scenario Scenario Scenario Scenario _ 
a b a-b a b a-b a b a-b 
A 46.3 71.1 -24.8 82.0 87.6 —5.6 57.6 70.1 -12.5 
Al 69.5 85.3 -15.9 95.5 91.2 4.3 55.9 80.6 -24.7 
A2 68.8 93.8 -25.0 49.3 92.0 —42.7 67.0 87.7 -20.7 
A3 19.3 49.8 -30.5 65.7 95.0 -29.3 47.1 63.0 -15.9 
A4 36.3 46.4 —10.1 61.1 89.4 -28.3 48.6 57.6 -9.0 
A5 B2 58.4 -26.3 9319 92.0 r9 48.3 57.4 -9.1 


Notes: A refers to arid and semi-arid regions of Asia; Al1—A5 are the Caspian Coastal Plain, the Indus River Plain, eastern Mongolian 


Plateau, the Ordos Plateau and eastern Kazakhskiy Melkosopochnik, respectively; scenario b (or a) refers to AOD (aerosol optical depth) 
is (or not) considered as the control variables in partial correlation between PAR and NDVI; and scenario a-b is the percentage 
difference between scenario a and b. Abbreviations are the same as in Tables 2—4. 


with significantly positive correlation of all regions was larger in scenario b than in scenario a, 
indicating that when excluding the influence of AOD, partial correlation between PAR and NDVI 
had increased in arid and semi-arid regions of Asia. By comparison, at the middle growing stage, 
there was little change in area percentage between two scenarios. 

The evolution of partial correlation between PAR and NDVI at different growing stages (Figs. 
5 and S1) showed that when excluding the influence of AOD (scenario b), area percentage with 
significantly positive correlation was found in the Caspian Coastal Plain and Turan Plain at the 
early growing stage, in eastern Kazakhskiy Melkosopochnik and southeastern Mongolian 
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Fig. 5 Spatial pattern of partial correlation between NDVI (normalized difference vegetation index) at growing 
stages and PAR (photosynthetically active radiation) when considering AOD (aerosol optical depth) as the control 
variable (right panel) or not (left panel) in arid and semi-arid regions of Asia during three selected 16-d periods of 
growing stages: 97-112 d (a and b), 193-208 d (c and d), and 209-224 d (e and f). Areas with insignificant 
(P>0.05) correlations and those with missing AOD values are not shown. 


JOURNAL OF ARID LAND 


Plateau at the middle growing stage, and southern Mongolian Plateau at the late growing stage. 
This phenomenon might be due to the obvious different responses of vegetation with different 
canopy structures on PAR at three growing stages. 


3.3.2 Effect of PRE on vegetation activity 


Partial correlation between PRE and NDVI in scenario b (i.e., excluding the influence of both 
AOD and other three climatic factors) showed that area percentage with significantly positive 
correlation was larger than that in scenario a (i.e., excluding the influence of other three climatic 
factors), especially at the middle growing stage (Table 2), implying the adverse impact of AOD 
on NDVI. Although five typical regions (A1—-A5) showed similar changes, there were much 
notable negative effects of atmospheric aerosols on vegetation growth over the Ordos Plateau 
(A4) at the early growing stage and over eastern Mongolian Plateau (A3) at the late growing 
stage, which is probably because atmospheric aerosols inhibited vegetation growth by 
suppressing precipitation via ‘burning’ the cloud droplets (Table 2; Figs. 6 and S2). Among five 
typical regions, there were still spatial differences in the response of vegetation activity at 
growing stages on precipitation when the effect of atmospheric aerosols was eliminated (i.e., 
changes from scenario a to scenario b). For instance, at the early growing stage, there was an 
evident increase in area percentage with a significant correlation (PRE vs. NDVI) over the Ordos 
Plateau when influence of AOD was excluded, whereas there was little change between two 
scenarios over eastern Mongolian Plateau. This difference at the early growing stage was possibly 
because of their different ways of water supply to vegetation growth, i.e., plants over the Ordos 
Plateau were strongly constrained by precipitation, and atmospheric aerosols could influence 
local vegetation by altering precipitation; while area percentage on eastern Mongolian Plateau 
was more affected by spring snow-melt water, thus leading the role of precipitation in vegetation 
activity less important so that the influence of AOD was diminished. By comparison, from the 
middle to late growing stage, there was an evident increase in areas with significant correlations 
over eastern Mongolian Plateau rather than the Ordos Plateau from scenarios a to b. During these 
stages, precipitation became the dominant factor limiting vegetation growth over eastern 
Mongolian Plateau due to its continental arid climate, whereas the Ordos Plateau was under the 
control of summer monsoon with relatively ample precipitation, making local vegetation less 
vulnerable to precipitation variations so that the influence of AOD on NDVI was relatively less 
important. 


Table 2 Area percentage (%) with significantly positive correlations from two scenarios of partial correlation 
between PRE (total precipitation) and NDVI (normalized difference vegetation index) in different regions at 
growing stages 


Early growing stage Middle growing stage Late growing stage 
Region Scenario Scenario Scenario Scenario Scenario Scenario Scenario Scenario Scenario 
a b a-b a b a-b a b a—b 
A 62.5 91.9 -29.4 50.9 86.3 -35.4 60.3 82.3 —22.0 
Al 50.3 87.8 -37.5 36.2 80.6 44.5 87.3 92.1 4.8 
A2 65.9 70.4 —4.5 0.7 SLE -90.5 40.6 58.8 -18.2 
A3 63.7 93.6 -29.9 50.3 90.3 —40.0 40.2 85.2 —45.0 
A4 26.7 76.7 —50.0 65.0 95.6 -30.6 76.6 92.6 —16.0 
A5 82.6 92.6 —10.0 373 81.7 —44.4 40.6 61.1 —20.6 


3.3.3 Effect of SSR on vegetation activity 

Partial correlation between SSR and NDVI showed that area percentage of the entire study area 
with significantly positive correlation was larger in scenario b than in scenario a, particularly at 
the early and middle growing stages (Table 3), indicating that the relationship between SSR and 
NDVI was closer without the influence of AOD. This was basically true for five typical regions, 
with the only exception of the Indus River Plain (A2). When AOD was used as the control 
variable (scenario b), NDVI was significantly and positively correlated with SSR over eastern 
Turan Plain and eastern Kazakhskiy Melkosopochnik at the early growing stage, and eastern and 
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central Mongolian Plateau at the late growing stage (Figs. 7 and S3). However, area percentage 
with a significantly positive correlation was smaller over the Indus River Plain (A2) and the 
Ordos Plateau (A3) at the early growing stage. 
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Fig. 6 Spatial pattern of partial correlation between NDVI (normalized difference vegetation index) at growing 
stages and PRE (total precipitation) when considering AOD (aerosol optical depth) as the control variable (right 
panel) or not (left panel) in arid and semi-arid regions of Asia during three selected 16-d periods of growing 
stages: 97-112 d (a and b), 193-208 d (c and d), and 209-224 d (e and f). Areas with insignificant (P>0.05) 


correlations and those with missing AOD values are not shown. 


Table 3 Area percentage (%) of significantly positive correlations from two scenarios of partial correlation 
between SSR (surface solar radiation) and NDVI (normalized difference vegetation index) in different regions at 
growing stages 


Early growing stage Middle growing stage Late growing stage 
Region Scenario Scenario Scenario Scenario Scenario Scenario Scenario Scenario Scenario 
a b a-b a b a—b a b a—b 
A 50.2 76.1 —26.0 17.9 44.9 -27.0 43.7 Se —10.0 
A1 48.1 82.5 -34.4 4.7 42.6 -37.9 43.6 48.5 —4.9 
A2 35.3 22.4 13.0 51.3 85.4 —34.1 36.4 57.4 -21.0 
A3 70.3 79.0 -8.7 222 85.5 —63.3 14.3 53.6 -39.3 
A4 41.6 54.4 -12.8 35.6 84.3 —48.7 60.8 79.4 -18.7 
A5 53.4 70.0 —16.6 5.6 23.4 -17.8 46.1 46.8 -0.7 


3.3.4 Effect of TEM on vegetation activity 


When the influence of aerosols was eliminated in partial correlation, area percentage with a 
significantly positive correlation between TEM and NDVI generally increased (i.e., from scenario 
a to b), which was more pronounced at the middle growing stage than at the other growing stages 
(Table 4). A similar change pattern was also found in five typical regions (Table 4). When AOD 
was used as the control variable that was eliminated in partial correlation, there were significantly 
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positive correlations over the Turan Plain, eastern Kazakhskiy Melkosopochnik and the Caspian 
Coastal Plain at the early growing stage; over northern Kazakhskiy Melkosopochnik and the 
Indus River Plain at the late growing stage (Figs. 8 and S4). 
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Fig. 7 Spatial pattern of partial correlation between NDVI (normalized difference vegetation index) and SSR 
(surface solar radiation) when considering AOD (aerosol optical depth) as the control variable (right panel) or not 
(left panel) in arid and semi-arid regions of Asia during three selected 16-d periods of growing stages: 97-112 d 
(a and b), 193-208 d (c and d), and 209-224 d (e and f). Areas with insignificant (P>0.05) correlations and those 
with missing AOD values are not shown. 


Table 4 Area percentage (%) of significantly positive correlations from two scenarios of partial correlation 
between TEM (surface air temperature) and NDVI (normalized difference vegetation index) in different regions at 
growing stages 


Early growing stage Middle growing stage Late growing stage 
Region Scenario Scenario Scenario Scenario Scenario Scenario Scenario Scenario Scenario 
a b a-b a b a-b a b a-b 
A 71.8 72:9 —2.1 19.0 40.3 —21.2 B83% 54.9 -19.3 
A1 49.9 80.5 -30.6 4.3 40.0 -35.6 SIL 64.6 -12.9 
A2 0.0 13.1 -13.1 19.8 49.2 -29.4 0.0 24.9 -24.9 
A3 75.6 91.9 -16.4 16.7 75.9 59.3 29.8 65.0 -35.1 
A4 79.0 75.3 3.7 10.4 69.7 59.3 41 30.2 —26.1 
AS 80.5 84.1 -3.6 4.0 14.4 -10.5 68.8 TS 4.9 


4 Discussion 


Atmospheric aerosols affect terrestrial vegetation activity indirectly via changing climate systems 
and ecosystems, including enhancing the scattering and absorption of solar radiation, changing 
atmospheric temperature, altering the probability and quantity of precipitation and providing 
nutrients to some extent (Ramanathan et al., 2001; Steiner et al., 2005; Wild, 2009). In addition, 
atmospheric aerosols may influence vegetation production and activity through the modulation of 
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Fig. 8 Spatial patterns of partial correlation between NDVI (normalized difference vegetation index) at growing 
stages and TEM (surface air temperature) when considering AOD (aerosol optical depth) as the control variable 
(right panel) or not (left panel) in arid and semi-arid regions of Asia of three selected 16-d periods at growing 
stages: 97-112 d (a and b), 193-208 d (c and d), and 209-224 d (e and f). Areas with insignificant (P>0.05) 
correlations and those with missing AOD values are not shown. 


their photosynthesis, respiration, metabolism and water utilization (Tie et al., 2016). In the 
present study, our results showed that atmospheric aerosols had clear effects on the vegetation 
activity by influencing SSR, PAR, TEM and PRE in arid and semi-arid regions of Asia. 


4.1 Relationship between major climatic factors and vegetation activity when excluding 
the influence of aerosols 


Over the arid and semi-arid regions of Asia, a large proportion of areas exhibited significantly 
positive correlations between major climatic factors and NDVI when the influence of AOD was 
excluded. Firstly, area percentage with a significantly positive correlation between vegetation 
activity and PAR became higher when the influence of AOD was excluded, which could be 
attributable to the decreases in total PAR that could have facilitated vegetation growth. There 
were two plausible mechanisms behind the decreased PAR: (1) One explanation was that in the 
majority of study area, the regions with sparsely covered vegetation had emitted relatively higher 
amount of dust aerosols and resultant decreasing PAR, which inhibited vegetation growth. And (2) 
the other explanation was that total PAR was still decreasing in areas with dense vegetation (such 
as northeastern Inner Mongolia) due to the high concentration of aerosols from their source areas 
and consequently the decreasing PAR reduced the areas in which PAR facilitates vegetation 
growth (Xie et al., 2010). 

Secondly, area percentage with a significantly positive correlation between SSR and NDVI 
became smaller when the influence of AOD was considered, which was possibly because of the 
negative forcing of higher aerosol emissions on SSR in arid and semi-arid regions of Asia with 
sparse vegetation coverage (Lau et al., 2006; Zhang et al., 2009; Huang et al., 2014). Generally, 
the surface cooling effects of aerosols (Du et al., 2014; Huang et al., 2014) would delay the onset 
of vegetation phenology in spring with reduced material accumulation of vegetation, and would 
trigger a higher risk of frost with restrained vegetation growth and decreased nutrient availability 
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at the entire growing stage (Bokhorst et al., 2008; Bennie et al., 2010; He et al., 2015). 

Thirdly, area percentage with a significantly positive correlation between PRE and NDVI was 
larger when eliminating the effects of AOD, possibly because aerosols may suppress precipitation 
by altering the cloud particle size and cloud temperature (Yin and Chen, 2007; Rosenfeld et al., 
2014; Huang et al., 2016), thus restraining vegetation activity in arid and semi-arid regions 
(Schimel et al., 2001; Fang et al., 2004; Chuai et al., 2013). Moreover, enhanced aerosol 
emissions from potential source areas were generally accompanied by reduced TEM and 
suppressed precipitation (Piao et al., 2006; Gessner et al., 2013), which would be unfavorable for 
vegetation activity. It should be noted that the Inner Mongolian Plateau did not show significantly 
positive correlation between NDVI and major climatic factors when AOD was used as a control 
variable, which was possibly because vegetation of this region might be more attributable to 
anthropogenic activities than to atmospheric aerosols. For instance, the afforestation program has 
been proven to promote vegetation restoration in recent years in China (He et al., 2015; Liu et al., 
2015). 


4.2 Spatiotemporal variation in the relationship between major climatic factors and 
vegetation activity 


There were obvious spatial and temporal differences in responses of vegetation to climatic factors 
under the influence of aerosols. For instance, at the early growing stage, there were significantly 
positive correlations between NDVI and climatic factors including PAR, SSR and TEM over the 
Caspian Coastal Plain and the Turan Plain when the influence of aerosols was excluded, which 
were all correlated with the higher concentration of spring aerosols that had substantially reduced 
PAR, SSR and TEM via altering the incident (diffuse and direct) solar radiation (Cheng et al., 
2015), leading to decreased contributions of these factors to vegetation growth (Mohammat et al., 
2013; Huang et al., 2014; He et al., 2015; Bi et al., 2017). 

At the middle growing stage, area percentages with significantly positive correlations of NDVI 
with PAR, SSR and TEM were only found in central and western Kazakhskiy Melkosopochnik 
and southeastern Mongolian Plateau when the influence of aerosols was not considered, which 
could be attributed to the following reasons: firstly, photosynthesis of vegetation over these 
regions with a simple canopy structure would easily reach its saturation points in summer, even 
though dramatic reduction of aerosol concentrations in these regions could increase both PAR and 
SSR (Niyogi et al., 2004; Strada et al., 2015); secondly, the increased surface temperature resulted 
from descending aerosol concentrations was not any more the key climatic factor influencing 
vegetation growth in hot summers (Gessner et al., 2013; Mohammat et al., 2013). However, there 
was still a significantly positive correlation between PRE and NDVI with an extra exclusion the 
influence of aerosols in most regions because dust aerosols (the main type of atmospheric 
aerosols) remained an important factor in suppressing PRE in summer that dust aerosols embed in 
the cloud layer would reduce cloud cover (especially the lower cloud) and water vapor content, 
and would increase the evaporation of cloud droplets, thus eventually constraining vegetation 
growth (Huang et al., 2006; Huang et al., 2014; Li et al., 2018). 

At the late growing stage, area percentages with significantly positive correlations between 
NDVI and PAR as well as SSR evidently increased over the Indus River Plain when excluding the 
influence of aerosols, possibly because the aerosols emitted from this region were primarily 
generated by anthropogenic activities, which would absorb SSR and PAR more easily than other 
types of aerosols (Dey et al., 2004). In addition, area percentage with a significantly positive 
correlation between TEM and NDVI also increased on the Indus River Plain probably due to the 
enhanced soil evapotranspiration driven by increased temperature, which was associated with an 
increased proportion of anthropogenic activities (Latha et al., 2018). By comparison, a 
significantly positive correlation between PRE and NDVI was observed in most of study area 
with the only exception of the Indus River Plain, possibly due to the suppressed precipitation that 
excessive inorganic (dust) aerosol nuclei could reduce the volume of droplets (Qian et al., 2009; 
Li et al., 2011; Rosenfeld et al., 2014). 


5 Conclusions 


From 2005 to 2015, there were decreasing trends in NDVI over the regions with sparse vegetation 
coverage, such as the Caspian Coastal Plain, which had a maximum decreasing rate of >1x10°7/a. 
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However, AOD tended to increase as NDVI decreased in most regions. When there was an extra 
exclusion of the influence of aerosols, the responses of vegetation to climatic factors showed clear 
spatiotemporal differences. Generally, at different growing stages, aerosols could influence 
vegetation activity by changing different climatic factors. At the early growing stage, the higher 
concentration of aerosols would reduce PAR, SSR and TEM, thereby constraining vegetation 
activity, which was obvious over northern Kazakhskiy Melkosopochnik. At the middle growing 
stage, decreased concentration of aerosols had little effects on PAR, SSR and TEM, but greatly 
suppressed precipitation and eventually regulated vegetation activity. At the late growing stage, 
the potential effects of aerosols on vegetation activity were evident in areas with high 
anthropogenic activities, such as in the Indus River Plain. Different responses of vegetation 
activity to aerosols and their related climatic factors would be due to differences in the vegetation 
canopy structures, their light utility efficiency, as well as regional environmental and climatic 
conditions. Our studies on the impact of aerosols on vegetation activity in arid and semi-arid 
regions of Asia would not only clarify the ecological effects of the aerosols in the source areas, 
but also explain the roles of aerosols in climate systems and ecosystems. 


Acknowledgements 


This work was supported by the National Key Research and Development Program of China (2016YFA0601900), 
the Key Frontier Program of the Chinese Academy of Sciences (QYZDJ-SSW-DQC043), the National Natural 
Science Foundation of China (41771012) and the Applied and Basic Research Program from Tangshan Science 
and Technology Bureau, China (20130202b). 


References 


Bennie J J, Wiltshire A J, Joyce A N, et al. 2010. Characterizing inter-annual variation in the spatial pattern of thermal 
microclimate in a UK upland using a combined empirical—physical model. Agricultural and Forest Meteorology, 150(1): 
12-19. 

Bi J R, Huang J P, Shi J S, et al. 2017. Measurement of scattering and absorption properties of dust aerosol in a Gobi farmland 
region of northwestern China—a potential anthropogenic influence. Atmospheric Chemistry and Physics, 17: 7775-7792. 

Bokhorst S, Bjerke J W, Bowles, F W, et al. 2008. Impacts of extreme winter warming in the sub-Arctic: growing season 
responses of dwarf shrub heathland. Global Change Biology, 14(11): 2603-2612. 

Chameides W L, Yu H, Liu S C, et al. 1999. Case study of the effects of atmospheric aerosols and regional haze on agriculture: 
An opportunity to enhance crop yields in China through emission controls. Proceedings of the National Academy of 
Sciences of the United States of America, 96(24): 13626-13633. 

Cheng S J, Bohrer G, Steiner A L, et al. 2015. Variations in the influence of diffuse light on gross primary productivity in 
temperate ecosystems. Agricultural and Forest Meteorology, 201: 98-110. 

Chuai X W, Huang X J, Wang W J, et al. 2013. NDVI, temperature and precipitation changes and their relationships with 
different vegetation types during 1998-2007 in Inner Mongolia, China. International Journal of Climatology, 33(7): 
1696-1706. 

Dey S, Tripathi S N, Singh R P. 2004. Influence of dust storms on the aerosol optical properties over the Indo—Gangetic Basin. 
Journal of Geophysical Research, 109: D20211. 

Dong B, Wilcox L J, Highwood E J, et al. 2019. Impacts of recent decadal changes in Asian aerosols on the East Asian summer 
monsoon: roles of aerosol-radiation and aerosol-cloud interactions. Climate Dynamics, 53(5): 3235-3256. 

Du J, He Z B, Yang J J, et al. 2014. Detecting the effects of climate change on canopy phenology in coniferous forests in 
semi-arid mountain regions of China. International Journal of Remote Sensing, 35(17): 6490-6507. 

Dubovik O, King M D. 2000. A flexible inversion algorithm for retrieval of aerosol optical properties from sun and sky 
radiance measurements. Journal of Geophysical Research, 105(D16): 20673-20696. 

Engström A, Ekman A M L. 2010. Impact of meteorological factors on the correlation between aerosol optical depth and cloud 
fraction. Geophysical Research Letters, 37(18): 1480-1493. 

Ezhova E, Ylivinkka I, Kuusk J, et al. 2018. Direct effect of aerosols on solar radiation and gross primary production in boreal 
and hemiboreal forests. Atmospheric Chemistry and Physics, 18: 17863-17881. 

Fang J Y, Piao S L, He J S, et al. 2004. Increasing terrestrial vegetation activity in China, 1982—1999. Science China Series C, 
47(3): 229-240. 

Gessner U, Naeimi V, Klein I, et al. 2013. The relationship between precipitation anomalies and satellite-derived vegetation 
activity in Central Asia. Global Planetary Change, 110: 74-87. 

Ginoux P, Prospero J M, Gill T E, et al. 2012. Global-scale attribution of anthropogenic and natural dust sources and their 


JOURNAL OF ARID LAND 


emission rates based on MODIS deep blue aerosol products. Review of Geophysics, 50(3): RG3005. 

Gu L, Pallardy S G, Hosman K P, et al. 2016. Impacts of precipitation variability on plant species and community water stress 
in a temperate deciduous forest in the central US. Agricultural and Forest Meteorology, 217: 120—136. 

He B, Chen A F, Wang H L, et al. 2015. Dynamic response of satellite—derived vegetation growth to climate change in the Three 
North Shelter Forest Region in China. Remote Sensing, 7(8): 9998-10016. 

Hsu N C, Gautam R, Sayer A M, et al. 2017. Global and regional trends of aerosol optical depth over land and ocean using 
SeaWiFS measurements from 1997 to 2010. Atmospheric Chemistry and Physics, 17: 8037-8053. 

Huang J P, Minnis B, Lin T, et al. 2006. Possible influences of Asian dust aerosols on cloud properties and radiative forcing 
observed from MODIS and CERES. Geophysical Research Letters, 33(6): L06824. 

Huang J P, Fu Q, Su J, et al. 2009. Taklimakan dust aerosol radiative heating derived from CALIPSO observations using the 
Fu—Liou radiation model with CERES constraints. Atmospheric Chemistry and Physics, 9: 4011—4021. 

Huang J P, Wang T H, Wang W C, et al. 2014. Climate effects of dust aerosols over East Asian arid and semiarid regions. 
Journal of Geophysical Research, 119: 11398-11416. 

Huang J P, Li Y, Fu C, et al. 2017. Dryland climate change: Recent progress and challenges. Review of Geophysics, 55(3): 
719-778. 

Hui W J, Cook B I, Ravi S, et al. 2008. Dust—rainfall feedbacks in the West African Sahel. Water Resources Research, 44(5): 
W05202. 

IPCC. 2013. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment 
Report of the Intergovernmental Panel on Climate Change. New York: Cambridge University Press, 1-12. 

Justice C O, Vermote E, Townshend J R G, et al. 1998. The moderate resolution imaging spectroradiometer (MODIS): Land 
remote sensing for global change research. IEEE Transactions on Geoscience and Remote Sensing, 36(4): 1228-1249. 

Kim Y, Kimball J S, Didan K, et al. 2014. Response of vegetation growth and productivity to spring climate indicators in the 
conterminous United States derived from satellite remote sensing data fusion. Agricultural and Forest Meteorology, 194: 
132-143. 

Latha R, Vinayak B, Murthy B S. 2018. Response of heterogeneous vegetation to aerosol radiative forcing over a northeast 
Indian station. Journal of Environmental Management, 206: 1224-1232. 

Lau K M, Kim M K, Kim K M, 2006. Asian summer monsoon anomalies induced by aerosol direct forcing: The role of the 
Tibetan Plateau. Climate Dynamics, 26: 855-864. 

Levy R C, Mattoo S, Munchak L A, et al. 2013. The collection 6 MODIS aerosol products over land and ocean. Atmospheric 
Measurement Techniques, 6: 2989-3034. 

Li X Q, Ting M F, Lee D E. 2018. Fast adjustments of the Asian summer monsoon to anthropogenic aerosols. Geophyscial 
Research Letters, 45(2): 1001- 1010. 

Li Z Q, Niu F, Fan J W, et al. 2011. Long-term impacts of aerosols on the vertical development of clouds and precipitation. 
Natural Geoscience, 4: 888-894. 

Liu X F, Zhu X F, Li S S, et al. 2015. Changes in growing season vegetation and their associated driving forces in China during 
2001-2012. Remote Sensing, 7(11): 15517-15535. 

Mi W, Li Z Q, Xia X G, et al. 2007. Evaluation of the Moderate Resolution Imaging Spectroradiometer aerosol products at two 

aerosol robotic network stations in China. Journal of Geophysical Research, 112: D22S08. 

Mike B. 2002. The Influence of Aerosols on Plant Growth. Honolulu: Georgia Institute of Technology, 5-10. 

Min Q L, Li R, Lin B, et al. 2009. Evidence of mineral dust altering cloud microphysics and precipitation. Atmospheric 

Chemistry and Physics, 9(9): 3223-3231. 

Mohammat A, Wang X H, Xu X T, et al. 2013. Drought and spring cooling induced recent decrease in vegetation growth in 

Inner Asia. Agricultural and Forest Meteorology, 178-179: 21-30. 

Nair U S, Mcnider R T, Patadia F, et al. 2011. Sensitivity of nocturnal boundary layer temperature to tropospheric aerosol 

surface radiative forcing under clear-sky conditions. Journal of Geophysical Research, 116(D2): D02205. 

Nemani R R, Keeling C D, Hashimoto H, et al. 2003. Climate—driven increases in global terrestrial net primary production from 

1982 to 1999. Science, 300(5625): 1560-1563. 

Niyogi D, Chang H I, Saxena V K, et al. 2004. Direct observations of the effects of aerosol loading on net ecosystem CO, 


exchanges over different landscapes. Geophysical Research Letters, 31(20): L20506. 

Park H S, Sohn B J. 2010. Recent trends in changes of vegetation over East Asia coupled with temperature and rainfall 
variations. Journal of Geophysical Research, 115: D14101. 

Pauli H, Gottfried M, Dullinger S, et al. 2012. Recent plant diversity changes on Europe's mountain summits. Science, 
336(6079): 353-355. 

Peng S, Chen A, Xu L, et al. 2011. Recent change of vegetation growth trend in China. Environmental Research Letters, 6(4): 
044027. 

Piao S L, Mohammat A, Fang J Y, et al. 2006. NDVI-based increase in growth of temperate grasslands and its responses to 
climate changes in China. Global Environmental Change, 16: 340-348. 


JIAO Linlin et al.: Potential responses of vegetation to atmospheric aerosols in arid and semi-arid regions... 


Qian Y, Gong D Y, Fan J W, et al. 2009. Heavy pollution suppresses light rain in China: Observations and modeling. Journal of 
Geophysical Research, 114: DOOKO2. 

Ramanathan V, Crutzen P, Kiehl J, et al. 2001. Aerosols, climate, and the hydrological cycle. Science, 294(5549): 2119-2124. 

Roderick M L, Farquhar G D, Berry S L, et al. 2001. On the direct effect of clouds and atmospheric particles on the 
productivity and structure of vegetation. Oecologia, 129(1): 21-30. 

Rosenfeld D, Rudich Y, Lahav R, et al. 2001. Desert dust suppressing precipitation: A possible desertification feedback loop. 
Proceedings of the National Academy of Sciences of the United States of America, 98(11): 5975-5980. 

Rosenfeld D, Wood R, Donner L J, et al. 2013. Aerosol cloud—mediated radiative forcing: Highly uncertain and opposite effects 
from shallow and deep clouds. New York: Springer Science+Business Media, 105-149. 

Rosenfeld D, Andreae M O, Asmi A, et al. 2014. Global observations of aerosol-cloud-precipitation-climate interactions. 
Reviews of Geophysics, 52(4): 750-808. 

Running S W, Justice C O, Salomonson V, et al. 1994. Terrestrial remote sensing science and algorithms planned for 
EOS/MODIS. International Journal of Remote Sensing, 15(17): 3587-3620. 

Schimel D S, House J I, Hibbard K A, et al. 2001. Recent patterns and mechanisms of carbon exchange by terrestrial 
ecosystems. Nature, 414: 169-172. 

Shao Y P, Wyrwoll K H, Chappell A, et al. 2011. Dust cycle: An emerging core theme in Earth system science. Aeolian 
Research, 2(4): 181-204. 

Solano R, Didan K, Jacobson A, et al 2010. MODIS Vegetation Index Users' Guide (MOD13 Series) Collection 5. Tucsn: 
University of Arizona, 35—42. 

Simmons A, Uppala S, Dee D, et al. 2006. ERA-interim: New ECMWF reanalysis products from 1989 onwards, ECMWF News, 
110: 25-35. 

Steiner A L, Chameides W L. 2005. Aerosol-induced thermal effects increase modelled terrestrial photosynthesis and 
transpiration. Tellus Series B: Chemical and Physical Meteorology, 57(5): 404—411. 

Strada S, Unger N, Yue X. 2015. Observed aerosol-induced radiative effect on plant productivity in the eastern United States. 
Atmospheric Environment, 122: 463—476. 

Tie X, Huang R J, Dai W, et al. 2016. Effect of heavy haze and aerosol pollution on rice and wheat productions in China. 
Scientific Reports, 6(1): 29612. 

Wang X, Wu J, Chen M, et al. 2018. Field evidences for the positive effects of aerosols on tree growth. Global Change Biology, 
24: 498—499. 

Wang X M, Ma W Y, Hua T, et al. 2017. Variation in vegetation greenness along China's land border. Science China Earth 
Sciences, 60: 1-8. 

Wild M. 2009. Global dimming and brightening: A review. Journal of Geophysical Research, 114: D00D16. 

Xie J, Huang J B, Wang G Y, et al. 2010. The effects of clouds and aerosols on net ecosystem CO, exchange over semi-arid 
Loess Plateau of Northwest China. Atmospheric Chemistry and Physics, 10: 8205-8218. 

Yan Z B, Lin Z H, Zhang H. 2014. The relationship between the East Asian subtropical westerly jet and summer precipitation 
over East Asia as simulated by the IAP AGCM4.0. Atmospheric Oceanic Science Letters, 7(6): 487—492. 

Yin Y, Chen L. 2007. The effects of heating by transported dust layers on cloud and precipitation: A numerical study. 
Atmospheric Chemistry and Physics, 7: 3497-3505. 

Yoon J, Burrows J P, Vountas M, et al. 2014. Changes in atmospheric aerosol loading retrieved from space-based measurements 
during the past decade. Atmospheric Chemistry and Physics, 14: 26001-26041. 

Yu H B, Chin M, West J J, et al. 2013. A multi-model assessment of the influence of regional anthropogenic emission reductions 
on aerosol direct radiative forcing and the role of intercontinental transport. Journal of Geophysical Research, 118(2): 
700-720. 

Yue X, Unger N, Harper K, et al. 2017. Ozone and haze pollution weakens net primary productivity in China. Atmospheric 
Chemistry and Physics, 17(9): 6073—6089. 

Zhang D F, Zakey A S, Gao X J, et al. 2009. Simulation of dust aerosol and its regional feedbacks over East Asia using a 
regional climate model. Atmospheric Chemistry and Physics, 9: 1095—1110. 

Zhang R, Ouyang Z T, Xie X, et al. 2016. Impact of climate change on vegetation growth in arid Northwest of China from 1982 
to 2011. Remote Sensing, 8(5): 364. 

Zhang X Q, Lei Y C, Pang Y, et al. 2014. Tree mortality in response to climate change induced drought across Beijing, China. 
Climatic Change, 124(1): 179-190. 

Zhang X Y, Arimoto R, An Z S. 1997. Dust emission from Chinese desert sources linked to variations in atmospheric circulation. 
Journal of Geophysical Research, 102: 28—41. 

Zhou L M, Tucker C J, Kaufmann R K, et al. 2001. Variations in northern vegetation activity inferred from satellite data of 
vegetation index during 1981 to 1999. Journal of Geophysical Research, 106: 20069-20083. 


202106.00010v1 


chinaXiv 


ChinaXivA ERAT 


JOURNAL OF ARID LAND 


Appendix 


40°E 60°E 80°E 100°E 120°E 40°E 60°E 80°E 100°E 120°E 


40°E 60°E 80°E 100°E 120°E 0°E 60°E 80°E 100°E 120°E 
Approved No. GS(2021)1992 
Number of significantly positive correlation 


0 1 2 3 4 5 6 


Fig. S1 Spatial pattern in the number of 16-d period with a significantly positive correlation between NDVI 
(normalized difference vegetation index) and PAR (photosynthetically active radiation) when considering AOD 
(aerosol optical depth) as the control variable (right panel) or not (left panel) over the arid and semi-arid regions 
of Asia during growing stages (each having four or six 16-d periods). Areas with insignificant (P>0.05) 
correlations and those with missing AOD values are not shown. Abbreviations are the same as in Figures S2-S4. 
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Fig. S2 Spatial pattern in the number of 16-d period with a significantly positive correlation between NDVI and 
PRE (total precipitation) when considering AOD as the control variable (right panel) or not (left panel) over the 
arid and semi-arid regions of Asia during growing stages (each having four or six 16-d periods). Areas with 
insignificant (P>0.05) correlations and those with missing AOD values are not shown. 
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Fig. S3 Spatial pattern in the number of 16-d period with a significantly positive correlation between NDVI and 
SSR (surface solar radiation) when considering AOD as the control variable (right panel) or not (left panel) over 
the arid and semi-arid regions of Asia during growing stages (each having four or six 16-d periods). Areas with 
insignificant (P>0.05) correlations and those with missing AOD values are not shown. 
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Fig. S4 Spatial pattern in the number of 16-d period with a significantly positive correlation between NDVI and 
TEM (surface air temperature) when considering AOD as the control variable (right panel) or not (left panel) over 
the arid and semi-arid regions of Asia during growing stages (each having four or six 16-d periods). Areas with 
insignificant (P>0.05) correlations and those with missing AOD values are not shown. 


